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Abstract 



We study instanton effects in QCD at very high baryon density. In this 



O . regime instantons are suppressed by a large power of {Aqcd/ fJ-), where Aqcd 

is the QCD scale parameter and /i is the baryon chemical potential. Instan- 
tons are nevertheless important because they contribute to several physical 
^ . observables that vanish to all orders in perturbative QCD. We study, in partic- 

ular, the chiral condensate and its contribution tuq^ m{'ij)'ip) to the masses 
. of Goldstone bosons in the CFL phase of QCD with Nf = 3 flavors. We find 

that at densities /) ~ (5 — 10)pO) where po is the density of nuclear matter, 
the result is dominated by large instantons and subject to considerable uncer- 
tainties. We suggest that these uncertainties can be addressed using lattice 
calculations of the instanton density and the pseudoscalar diquark mass in 
QCD with two colors. We study the topological susceptibility and Witten- 
Veneziano type mass relations in both Nc = 2 and Nc = 3 QCD. 
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I. INTRODUCTION 



Strange quark matter at very high baryon density but low temperature is beheved to be 
in the color-flavor-locked (CFL) phase The CFL phase is characterized by diquark 

condensates which break both the global SU{3)l x SU{3)ii flavor symmetry and the local 
SU{3) color symmetry but preserve a global SU{3)v- In addition to that, diquark condensa- 
tion spontaneously breaks the exact U{l)v and approximate U{1)a symmetry of QCD. The 
low energy behavior of the CFL phase is governed by the corresponding Goldstone bosons, 
an octet of pseudoscalar mesons associated with chiral symmetry breaking and two singlets 
associated with the breaking of U{1)v,a- 

The CFL phase has many similarities with Nf = 3 QCD at low density but the 

mechanism of chiral symmetry breaking is quite different. In particular, if the baryon den- 
sity is very large, the dominant order parameter for chiral symmetry breaking is not the 
usual quark- ant i- quark condensate {ipip)-, but a four-fermion operator {{ipipY) ~ {ipip)'^. As 
a consequence, the masses squared of the pseudoscalar Goldstone bosons are quadratic, not 
linear, in the quark masses. The coefficient of proportionality can be computed in pertur- 
bative QCD g-Q. 

The CFL diquark condensate leaves a discrete axial {Z2)a symmetry unbroken. Since the 
quark- ant i- quark condensate is odd under this symmetry, {ipip) remains zero in weak coupling 
perturbation theory. Non-perturbative effects, instantons, break the {Z2)a symmetry and 
lead to a non-zero expectation value (ipip). If the density is large instantons are strongly 
suppressed and (ipip) is small. As the density decreases instantons become more important 
and the chiral condensate grows. Our main objective in the present work is to compute the 
linear term in the mass relation for the Goldstone bosons. The relative size of the linear and 
quadratic terms in the mass relation provides an estimate of the baryon density at which 
the crossover from the asymptotic regime, in which chiral symmetry breaking is dominated 
by the diquark condensate, to the low density regime, in which chiral symmetry breaking is 
governed by the quark-anti-quark condensate, occurs. 
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The size of the hnear term in the Goldstone boson mass relation also has important 
consequences for the phase structure of CFL matter in the presence of a non-zero strange 
quark mass and lepton chemical potentials p!3|- p!9|j . Bedaque and Schafer argued that for 
physical values of the quark masses and the baryon chemical potential CFL quark matter is 
likely to be kaon condensed [0. This conclusion was based on an analysis of the effective 
potential for the chiral order parameter. The low energy constants that appear in the 
effective potential were obtained in perturbation theory, but the conclusions are unchanged 
if the coefficients are estimated using dimensional analysis. This suggests that the results of 
the perturbative calculation are qualitatively correct even in the regime of strong coupling. 
However, the assumption that instanton effects are small for the densities of interest is 
crucial. 

The paper is organized as follows. In section H we compute the instanton contribution to 
Goldstone boson masses in the CFL phase. In section |T| we study corrections to this result 
due to the finite size of instantons. In section ^ we compute the pseudoscalar diquark mass 
in QCD with two colors. This observable can be studied with present lattice techniques. We 
also comment on Witten-Veneziano relations in QCD at high baryon density. Our results 
extend previous work in PD|,P|,|8|PT[ . 



II. INSTANTON CONTRIBUTION TO CFL CHIRAL THEORY 



In this section we consider the CFL phase of high density quark matter. For excitation 
energies smaller than the gap the only relevant degrees of freedom are the Goldstone modes 
associated with the breaking of chiral symmetry and baryon number. The interaction of the 
Goldstone modes is described by the effective Lagrangian 
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Here S = exp(i0'^A°/ /,r) is the chiral field, f-^ and / are the octet and singlet decay constants, 

M is a complex mass matrix. The chiral field and the mass matrix transform as S ^ LT,R^ 

and M — > LMR^ under chiral transformations (L, R) G SU{3)l x SU{3)[i. The axial U{1)a 

field is V = exp(i0) = exp(2i?7'/(\/6/)). As explained in ||T^ the covariant derivative VqS 

contains the effective chemical potentials Xl = {MM^)/{2^) and Xr = (MtM)/(2/i), 

(MM^\ (M^M\ 
V„E = a,E + . — E-ffi — . (2) 



The coefficients Bi of the quadratic mass terms were computed in [§-[T2| by matching the 
mass dependence of the vacuum energy computed in perturbative QCD and the CFL chiral 
theory. The result is P,p!2 



B, = -B, = ^, B, = 0, (3) 

where A is the gap in the quasi-particle spectrum. We shall determine the coefficient A by 
computing the shift in the vacuum energy linear in the quark mass. This shift is due to the 
first order instanton contribution to the vacuum energy in the background of a perturbatively 
generated diquark condensate, see the diagram shown in Fig. |l]a. In QCD with three fiavors, 
the instanton induced interaction between quarks is given by p3|- p5|] 

C = J n{p, fl)dp ^ i^ ^fif2h^9i92g^ (^^^(V'iJ,/iV'L,9i)(V'iJ,/2^L,92)(V'iJ,/3V'L,<;3) (4) 

-ii^R,h'ipL,9^){'^pR,f2(^^,u'^pL,92){'4'R,h(r^ui^L,g3) + {L ^ R)] . 



8{N, + 2) 

Here, p is the instanton size, p is the quark chemical potential, fi,gi are fiavor indices and 
o'fj.u = ^[ifj-ylu]- The instanton size distribution n{p, p) is given by 



n{p,p) = Cn f ^ ) p ^exp 



9{p) 



exp 



Nfp'p' , (5) 



_ 0.466 exp(-l. 679 ATJ 1.34^/ , 
- (iV.-l)!(iV.-2)! ' 

^""^ 61og(pA), h=^-^N,-lNf. (7) 



g\p) ^' 3^3 

At zero density, the p integral in equ. is divergent at large p. This is the well-known 
infrared problem of the semi-classical approximation in QCD. At large chemical potential, 



however, large instantons are suppressed and the typical instanton size isp~/i ^. 
To linear order in the quark mass one of the three zero modes is lifted. We obtain 

C = J n{p, fi)dp ^^fy^ef,f,f,eg,g,g.,Mf^g^ (^^^^(V'iJ,/iV'L,3i)(V'K,/2V'L,92) (8) 

i'ipR,hCr^^lljL,g,){lpR,f2^l„y'^L,g2) + {M ^ M\ L ^ R) 



We can now compute the expectation value of equ. (|[) in the CFL ground state P, pO| , pT 
Using the perturbative result for the diquark condensate in the CFL phase, 



(9) 
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we find the instanton contribution to the vacuum energy density 
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We note that for M = diag(mu, m^, ms) the instanton contribution to the vacuum en- 
ergy is indeed negative. Since the effective interaction involves both left and right-handed 
fermions the relative phase between the left and right-handed condensate in equ. is im- 
portant. Instantons favor the state with {ipL'^Pi) = —{4'r4'r) which is the parity even ground 
state. Equation ( p^ for the vacuum energy can be matched against the effective lagrangian 
equ. (p. We find 



A = C 



N- 



87r^r(6) 



W2n^ ( /i^ 
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where we have performed the integral over the instanton size p using the one-loop beta 
function. We note that A is related to the quark-anti-quark condensate, {ipip) = —2A. 



Equation ([TT|) agrees, up to a numerical factor and a power of g, with the estimate presented 
in ||2|J^. We can also determine the masses of Goldstone bosons. We take into account the 
instanton contribution equ. (0), the O(M^) term given in equ. and the O(M^) term 
given in equ. (Ef). To this order, the masses of the charged Goldstone bosons are given by 
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In the flavor symmetric limit m„ 
mass of the 1]' is 



2A 



4B 



J TV J TV 
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(12) 



= m the one-instanton contribution to the 
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-m. 



(13) 



If flavor symmetry is broken the rj' mixes with the T] and 7r° We note that at 

the one-instanton level, the mass of the rj' vanishes in the chiral limit. The r/'-mass in the 
chiral limit arises from the two-instanton diagrams shown in Figs. |l|c)-|l]e). These diagrams 
are hard to evaluate and we will not pursue this problem here. However, even without a 
calculation we can determine the dependence of the instanton generated potential on the 
QCD theta angle 6 and the U{1)a phase of the chiral field. Again restricting ourselves to 
the case of exact flavor symmetry we find 



V = -6mA cos{e + 0) - 12m^B cos(0) - 2C cos(20 + 30) , 



(14) 



where A is the coefficient of the one-instanton contribution given in equ. (|TlD , B = Bi = B2 
is the coefficient of the O(M^) term given in equ. (§), and C is the coefficient of the two- 
instanton term. The potential equ. (|l^ determines the topological susceptibility in the CFL 
phase. In the limit of very small quark masses we find 

2mA m{%l)ip) 



Xtop 



(15) 
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This result agrees with prediction of anomalous Ward identities at zero density. If we take 
the chemical potential to infinity while keeping the quark mass fixed then m'^B ^ mA ^ C 
and the topological susceptibility is given by xtop = GmA = —3m{ijjip). 

In the following, we will study the mass of the in more detail. The mass of the 
receives the largest instanton contribution and is of special interest in connection with 
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kaon condensation. Different contributions to tlie i^^-mass are sliow in Fig. |^. Tlie 0{m) 
and 0{m^) contributions are computed from equ. (0) taking into account only tlie terms 
proportional to A and B, respectively. The O(m^) term (m^ — m^)/(2/i) increases the mass 
of the , but decreases the mass of the Kq. The complete result for the mass is show in 
Fig. 1^. If the O(m^) term is bigger than the 0{m,im?) contribution then condensation 
takes place and the mass of the vanishes. We have used m„ = 4 MeV, = 7 MeV, 



150 MeV and the perturbative result for the gap Pq -p9 



A = 5127r^2-^/3(2/3)-^/'6[,/i(/-^exp 



(16) 



where b'^ is a constant which is determined by non-Fermi liquid effects [pO| , ^ that are not 
included in our calculation. To leading order in perturbation theory the pion decay constant 
is given by p|, p!0| , |32| , p3 

21 -81og(2) / \ 
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Qualitatively, the instanton contribution to the kaon mass scales as 



(17) 



rriK 
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7/2 



A. 



(18) 



A / VA, 

This shows that the result is strongly suppressed as /i — > oo. We also note, however, 
that the result is quite sensitive to the value of the scale parameter A. In practice, the 
power dependence on the scale parameter is canceled to some degree by the logarithmic 
dependence. This can be seen from the results shown in in Fig. 0. At /i = 500 MeV, the 



instanton contribution to mj^o calculated from equ. (|TT]) varies between 85 MeV and 120 
MeV if the scale parameter is varied between 180 MeV and 280 MeV. 

The dependence of ttiko on the scale parameter provides a naive estimate of the impor- 
tance of higher order corrections. In the present case this is probably an underestimate. This 
can be seen by studying the role of higher order corrections in the instanton size distribution. 
At two-loop order we have 



n//(p,/i) = Cjv (/5/(p))^^'p ^exp -f3ii{p) - Nfp^fi 



.2 ,.2 



(19) 



-61og(pA), 
b' 
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(20) 



(21) 



The instanton contribution to the kaon mass calculated with the two-loop instanton dis- 
tribution is also shown in Fig. ^ We observe that the results are significantly smaller as 
compared to the leading order estimate. We find mKo{inst) = (17 — 40) MeV compared 
to the leading order result mxo{inst) = (85 — 120) MeV. The large difference between the 
one and two-loop results is related to the fact that at moderate density perturbation theory 
predicts that the average instanton size is not small. Using equns. (|5|, p^) we get 



For /i = 500 MeV we find p = 0.55 fm, which is bigger than the standard estimate for 
the typical instanton size at zero density po = (0.3 — 0.4) fm. The situation is somewhat 
improved for the two-loop size distribution which gives p = 0.45 fm. The main difference 
between the one and two-loop size distributions is that at two-loop order the pre-exponent 
is evaluated at a scale given by the inverse instanton size p~^ rather than at the external 
scale /i. This difference is formally of higher order, but at moderate density it provides a 
significant suppression of large instantons. 

We can also study this problem in a different way. It is clear that at small density 
there has to be some non-perturbative effect that eliminates the contribution of large-size 
instantons. We can simulate this effect in terms of a non-perturbative screening factor 
exp(— p^m^^j.) in the instanton size distribution. If the screening mass is adjusted in such a 
way that the average instanton size at p = 500 MeV is equal to the phenomenological value 
at p = 0, p = 0.35 fm, then we find rriKoiinst) = (7 — 12) MeV. 

The sensitivity of the instanton contribution to the kaon mass to large-size instantons 
translates into a large uncertainty regarding the behavior of kaons at finite density. If large 
instantons with p ~ 0.5 fm play a role then the kaon mass is dominated by the instanton 
contribution and kaon condensation is unlikely. If the typical instanton size satisfies p < 




(22) 
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0.35 fm then the instanton contribution to the kaon mass is at most comparable to the 
perturbative contribution and kaon condensation is hkely. 



III. FINITE INSTANTON SIZE EFFECTS 

In the previous section we computed the instanton contribution to the vacuum energy 
using the effective instanton induced interaction equ. (|^). This effective interaction is derived 
under the assumption that the relevant momenta are smaller than the inverse instanton size, 
p,li <^ p^^. In our case the relevant momenta are on the order of the Fermi momentum 
IpI ~ ~ /i and the typical instanton size is given by p ~ /x^^, so pp ~ 1. This implies 
that finite size effects are potentially important. 

Instanton finite size effects can be determined from the fermion zero mode solution at 
finite baryon density. This solution was found in 0] and studied in |]35| , ^|36| . Using the 



exact zero mode solution amounts to the replacement 

i^fo.^ - ^i^fa^n. V'^"^ - (23) 

/vr /vr ■' 

in the effective interaction equ. (H). Here, the instanton form factors T'^ are given by 
2g,|3| 



n = ^- ■ crm. = liv* ■ <y~)(P ■ <y^% (24) 

with cr^ = (±ia, 1), p = (p,P4 + «/i), v^a = V^a(p, /^) and = [ipa{p, The instanton 

form factors are determined by the Fourier transform ipa{p, p^) of the instanton zero mode 
wave function. For completeness, we provide the result for (pa{p, P-) in appendix A. 

We first study how the instanton form factor modifies the loop integral which contains 
the mass insertion, see equ. (^. In this case we need the quark-anti-quark form factor 

Giip) = (P)-^I(P). (25) 

We find Gi{p) = {p + ipiYipaip^ p)Va{p^ p) / {'^t^Y ■ We can now perform the loop integral 
with the mass insertion. We find 



where we have used the normahzation condition equ. ( |A6| ) for the Fourier transform of 
the zero mode wave function. We observe that the instanton form factor does not modify 
equ. (1). 

The next step is the integration over the quark-quark propagators, equ. (p^Of) . The diquark 
loop involves the form factors 

FM = ^2^'^l^V)J'[{-vy,i (27) 

^2(P) = ^^e^'^(p)^^(-p)e,n.[(ao)(a -p)]!". (28) 

These two structures arise from the contraction of the instanton vertex with the diquark 
propagator 

S.(p) = i(l + 5-p)C7.;;j^^^^^^. (29) 

The two form factors -Fi,2(p) are given by 

F\{V) = -^^^{(P^ + /^^)V^(P'/^) ■ '^{-P,f^) + {p + if^) ■ V{P, f^){p - if^) ■ <^{-P,f^) 

- (p + if^) ■ ^{-P, f^){p - ifJ') ■viP^fJ')}, (30) 

If r ~* 

F2{p) = -^^-^[{p + ifi) ■ <^{p,fi) [{Pi - ifi)p ■ (p{-p, fi) -p^i{-p,fi) 

-{p- if^) ■ v{-Pi /^) [{P4 + ifj')p ■ 4>{p, /i) - m(p> ^A\ }• (31) 

The functions F±{p) = Fi{p) ± F2{p) are shown in Figs. H,^. We observe that the form 
factors are centered around the Fermi surface, P4 = 0, |p| = /i. This is an important 
observation. Instantons produce (t/'/jt/'l) -particle- hole pairs near the Fermi surface, and the 



chirality violating pair creation amplitude is therefore not suppressed by Pauli-blocking 
In Fig. ^ we study the dependence of F± [p) on p at fixed baryon chemical potential /i. We 
find that on the Fermi surface = 1 for all p whereas F_ ^ 1 only as p ^ 0. In practice 
we are interested in the limit /i — » oo with p ~ p^^. This limit is studied in Fig. ^. Again 
we observe that F+(p = /i,p4 = 0) = 1 for all p. 
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In the weak coupling limit the contraction of the diquark propagator equ. ( p9D with 
the instanton vertex only involves the form factor. Expanding around the Fermi 

surface and using 0) = 1 we can show that up to corrections of order 0{g^) the 

instanton form factor does not modify the weak coupling result equ. (|1T]) for A. In the same 
fashion we can also show that up to corrections of order 0{g^) there is no anti-gap 
contribution to A. 

We have also studied the effect of the instanton form factor at moderate density. The 
integral over the diquark propagator is given by 

f d^p f+(p)A(|p1) 

J (27r)>l + (|p1-/i)2 + A2- ^ ^ 

For F^{p) = 1 we recover the perturbative result equ. Results for the kaon mass 

with the effects of the instanton form factor included are shown in Fig. At a baryon 
chemical potential /i = 500 MeV the instanton form factor leads to a modest reduction of 
the instanton contribution to the kaon mass on the order of 40%. This reduction is due to 
the fact that the instanton form factor suppresses the contribution from large instantons. 
This effect, however, does not lead to a significant reduction of the average instanton size. 
In particular, the discrepancy between the one and two-loop results remains large. 



IV. QCD WITH TWO COLORS 

Given the large uncertainty in the instanton contribution to the Goldstone boson masses 
at moderate baryon density it would clearly be useful if the result could be checked on the 
lattice. Because of the sign problem direct studies of A^^^ = 3 QCD are still not feasible. On 
the other hand, QCD with Nc = 2 colors and Nf = 2 flavors does not suffer from a sign 



problem at non-zero baryon density and has been studied successfully on the lattice p8|-^ 

In QCD with two colors we also expect the formation of a diquark condensate at large 
baryon density ||3^. For Nf = 2 flavors the diquark condensate breaks the U{1)b of baryon 



number and the anomalous U{1)a symmetry, but not the SU (2)^ x SU {2)ji chiral symmetry. 
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The U{1)b Goldstone boson is exactly massless, but the U{1)a Goldstone boson acquires a 
mass from instantons. The corresponding contribution to the vacuum energy is shown in 
Fig. |I]f. Except for the mass insertion, the diagram is identical to the 0(m) contribution to 
the vacuum energy in the CFL phase of Nc = Nf = 3 QCD. 

The effective lagrangian for the pseudoscalar Goldstone boson is 



(33) 



We can think of as the relative phase between the left and right handed diquark condensate. 
Since baryon number is broken, the field has the quantum numbers of both pseudoscalar 
mesons and diquarks. This means that the mass of the 0-field governs the asymptotic 
behavior of both the rj' and pseudoscalar diquark correlation functions. The decay constant 
/ and Goldstone boson velocity v can be determined in perturbation theory. The result is 
identical to the one in QCD with Nc = 3 colors and Nf = 2 flavors |]5|,p!0[| 



47r2 



1 



(34) 



The potential V{(f)) is determined by instantons. The calculation of the instanton induced 
potential is completely analogous to the case N^ = 3, Nf = 2 so we can be brief here. 
We find V{(f)) = cos(0 — ff) where Q is the QCD theta angle and 



^2 = C^Stt^ 



4-^ i /^^ 



9 V2^l 
In A^c = 2 QCD the gap A is given by 

A = 5127rVo/i^"^exp 



A 



-2 



(35) 
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(36) 



Using these results we can determine the mass of the pseudoscalar Goldstone boson. We 
find 



P 



(37) 



This result looks like the Witten-Veneziano relation ||4l| , |42[| 
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= (38) 

Jr)' 

where A2 plays the role of the topological susceptibility^. We note, however, that the 
topological susceptibility vanishes in the chiral limit, Xtop = 0{rn?). Indeed, A2 is not the 
topological susceptibility but the density of instantons, A2 = {N/V). This is consistent with 
the idea that in a dilute system of instantons local fluctuations of the topological charge are 
Poissonian. This means that Xtop{V) = {Q"^) /V = {N/V) for any 4-volume V which is large 
compared to the average volume per instanton and small compared to the screening volume 

{N/V)-^ < y < r|„ (39) 

where r£) ~ m~,^ is the screening length. Because of the large mass of the rj' this window 
does not exist in QCD at zero density. In QCD at asymptotically high baryon density the 
decay constant scales as the baryon chemical potential, / ~ /i, the rj' mass is much smaller 
than {N/Vy/'^, and the window described by equ. ( pQ]) opens up pT |. 



Mass relations of the Witten-Veneziano type can also be derived in QCD with = ?> 
colors. The case of Nj = 2 flavors is exactly analogous to the Nc = Nf = 2 case except that 
the numerical value of the finite volume susceptibility XtopiV) is different |21|. In the Nf = 3 



CFL phase we also find that the r^'-mass in the chiral limit satisfies a finite-volume Witten- 
Veneziano relation, m'^,f^, = QxtopiY)- In this case, however, we have XtopiV) = 2{N/V) 
because instantons occur in pairs. 

We can also study how the relation between the 1]' mass and the topological susceptibility 
is modified if quark masses are included. In the regime in which the rj' mass is dominated by 
the linear mass term, equ. (|13]), we find rn^if^, = {4:/9)xtop{V). In this regime the topological 



^The extra factor 2Nf is related to our normalization of /. If we define /^/ in such a way that the 
canonical axial current Aq = V'7o75^V' is represented in the effective theory by Aq = fn'dQVj' + . . . 



then we find rnf^,f^, = 4^2- 
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charge is only partially screened, XtopiY c>o) = ^XtopiY)- If the quark mass is even larger 
and the quadratic mass term dominates the relation between the r/' mass and the topological 
susceptibility is lost. 

In Figs. |^,|^ we show numerical results for the rj' mass and the instanton density in 
QCD with Nc = Nf = 2. We again find that at moderate densities, jU ~ 500 MeV, the 
result is dominated by relatively large instantons p ~ 0.5 fm. The one-loop results for 
both the rj' mass and the instanton density are surprisingly large. Both are bigger than 
the phenomenological values m^/ = 945 MeV and (N/V) ^ lfm~^ in A^^c = 3 QCD at 
zero baryon density. This is incompatible with the idea that m^' and (N/V) decrease as a 
function of density, and that the dependence on Nc is weak. Lattice simulations of A^^c = 2 



also point to a very light pseudoscalar diquark P9,40|. We should note, however, that these 



simulations were performed with more than two continuum quark fiavors. 



V. SUMMARY 



We have computed the instanton contribution to the masses of Goldstone bosons in the 



CFL phase of QCD with three quark fiavors. Our main result is given in equns. (|Tl|) and 
(|T^). At very large baryon density this result is expected to be exact. 

At densities p ~ (5 — 10)po that are of interest in connection with the physics of neutron 
stars the perturbative result is very sensitive to the contribution of large instantons. For 
fi = 500 MeV leading order perturbation theory predicts p = 0.55 fm and mK{inst) = 
(85 — 120) MeV. If this result is reliable, then instantons would almost certainly prevent 
kaon condensation in the CFL phase at densities p < lOpo- The situation is different if the 
two- loop instanton size distribution is used. In this case we find p = 0.45 fm and mK{inst) = 
(17 — 40) MeV. The result is further reduced if a phenomenological screening mass and the 
instanton form factor is taken into account. In this case we find mxiinst) = (7 — 12) MeV. 

We suggest that these uncertainties can be addressed using lattice calculations of the 
instanton density and the pseudoscalar diquark mass in A'c = 2 QCD. The result of the 
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leading order calculation, equ. (|38|) suggests that the pseudoscalar diquark mass is much 
bigger than 2A. This prediction is only weakly dependent on the value of the scale parameter 
and the magnitude of the gap A. We also emphasize that the pseudoscalar diquark mass is 
related to the topological susceptibility in a finite volume. 
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APPENDIX A: FOURIER TRANSFORM OF FERMION ZERO MODES 



We repeat the results for the Fourier transform of the fermion zero mode wave function 
2^,36]. We follow the notation of |36|. The Fourier transform is given by 



The temporal and spatial components of are given by 

<^4(P4,P, /i) = ^1 - /i - ipi) [{2pi + ill) /i_ + i (p - /i - ipA) f2-] 



+ {p + 11 + ipi) [(2p4 + in) /i+ - i (p + /X + ipi) /s+j | 



<^(P4,P,/^) 



T^P Pi 

Ap 



{2p -p){p- p- ipi) fi- + (2p + /i) (p + /i + ipi) /: 



1+ 



1 



+ (^2(p - fi){p- fi- ipi) - - (/^ + m) [pI + {p~ ^)^] J /2- 

+ (2{p + n){p + fi + ip^) + ^ (/i + ^pa) [pI + (p + yu)'] j /2+ 

where p = p = p/p and we have introduced the functions 



(Al) 



(A2) 



(A3) 



/i± = — [h{z±)Ko{z±) - h{z±)K,{z±)] 



/2± = —h{z±)Ki{z^). 

4 



(A4) 
(A5) 



The argument of /i,2± is given by z± = \p\Jp'l + (p ± /i)^- The Fourier transform (f^{p, fx) 
satisfies the normalization condition 



d^p 



Vu{P, P)'^u{p, P) 



^b,(p, -/i)]v(p,/i) = ^. 



(A6) 
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FIGURES 




FIG. 1. Figure a) shows the 0(M) instanton contribution to the vacuum energy in the CFL 
phase. The six-fermion vertex is the effective 't Hooft vertex in the field of an instanton. Figure 
b) shows the O(M^) perturbative contribution to the vacuum energy. The four-fermion vertex 
corresponds to hard (p ~ pp) electric gluon exchange. Figures c),d) and e) show the two-instanton 
contribution to the vacuum energy in the chiral limit. Figure f) shows the instanton contribution 
to the vacuum energy in QCD with two colors and two flavors. 
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FIG. 2. Different contributions to the K mass at finite baryon density. Tfie two sets of dotted 
lines bounded by solid lines marked with triangles show the instanton contribution to the kaon 
mass computed for different values of the A parameter in the range A = (180 — 280) MeV. The 
curves marked with solid (open) triangles are computed with the one-loop (two-loop) instanton 
distribution. The solid lines marked with circles show the pcrturbative O(m^) contribution for 
two different values of the scale parameter. The solid line marked by a square shows the absolute 
magnitude of the effective chemical potential 0(/Xg^j) = O(m^) contribution. 
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11 [MeV] 



FIG. 3. This figure shows the sum of the different contributions to the K mass shown in 
Fig. |2[ For comparison we also show the gap 2A in the excitation spectrum for different values of 
the scale parameter in the range A = (180 — 280) MeV. 
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P [GeV] 

FIG. 4. Instanton form factors F± for quark-quark scattering. The form factors are shown 
as a function of p = \p\ for p4 = 0. The Fermi momentum was chosen to be pp = 0.5 GeV. The 
different curves correspond to instanton sizes in the range p = (0.15 — 0.90) fm. 



23 
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0.4 fm 
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- - n = 
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0.5 1 1.5 2 2.5 3 



p [GeV] 

FIG. 5. Same as Fig. |^ for different Fermi momenta pp = 0.5, 1.0, 2.0 GeV. The upper and 
lower curves show F- and F^, respectively. The solid square shows the value of -F+ on the Fermi 
surface. The instanton size was fixed at ppp = 1, corresponding to p = 0.4, 0.2, 0.1 fm. 
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finite size effects included 
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FIG. 6. Same as Fig. |^ but with instanton finite size effects (instanton form factors) included. 
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FIG. 7. Pseudoscalar Goldstone boson mass in A'^c = 2 QCD. The two sets of curves correspond 
to the one and two-loop instanton size distribution and different values of the scale parameter, see 
Fig. |2[ For comparison, we also show the energy gap 2A. 
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FIG. 8. Instanton density in Nc = 2 QCD. The two sets of curves correspond to the one 
and two-loop instanton size distribution and different values of the scale parameter in the range 
A = (180 - 280) MeV. 
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